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ABSTRACT 


Polycrystalline alloys Of iron in vanadium up to 15.32 
atomic percent were prepared principally by arc melting. In 
order to obtain the elastic constants and investigate the 
possibility of the existence of an elastically isotropic 
alloy ultrasonic measurements were made. These measurements 
determined the longitudinal and transverse wave velocities 
ise the 10-30 MEz range at room temperature. The resulting 
longitudinal modulus increased with alloying whereas the 
shear modunnae remained roughly constant although the data 
are somewhat scattered. These data are interpreted to show 
that the shear constants vary through alloying in a manner 
consistent with the results of other experiments and that 
the results are consistent with a s! da” as electronic 
configuration in the transition metals. The existence of 
an isotropic alloy is indicated for alloys with electron to 
atom ratios lower than vanadium rather than for the alloys 


measured. 
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I. ENTRODUCT TCH 


The decision to investigate the elastic properties of a 
vanadium-iron alloy was the result of a desire to have an 
elastically isotropic material such as tungsten, but without 
tungsten's high melting point, for use ina thin film 
experiment. It was felt that the isotropic property of such 
a material would also be advantageous in macroscopic poly- 
crystalline form. Isotropic grains would tend to deform 
uniformly in the presence of an overall stress irregardless 
of their orientations. This effect would reduce one source 
of motion at the grain boundaries resulting in a corresponding 
improved resistence to cyclic fatigue. 

Tt was known that the ratio of the elastic shear moduli, 
A = C/C', called the Zener anisotropy is less than unity for 
some of the bcc transition elements and greater than unity 
for others. Using a simple Vegards's Law model for the 
variation of A with alloying, various phase diagrams were 
investigated to find a pair of constituents where the re- 
quired relative concentrations for isotropy fell in a region 
of solid solubility of the system. The vanadium-iron system 
was found to have a solid solution bcc phase up to approxi- 
mately 40 at.% iron in vanadium. On the basis of Vegard's 
Law the concentration of iron for isotropy was estimated to 
be approximately 12 at.%. Another influencing factor in the 
choice of the V-Fe system was the recent availability of 


wery high pumity vanadium. 





The free atom electron configuration of vanadium is 
4s°3d° and of Leon, 4s*3a°. Thus the addi Elion of renee 
vanadium increases the electron to atom ratio of the alloy. 
According to Fuchs [Ref. 1] the resistance to shear in metals 
is a result of two interactions: (a) the Coulomb attraction 
between positive ions and negative conduction electrons, and 
(b) ene exchange interaction between ions, thought to be 
significant only between nearest and second nearest neighbors. 
This second interaction devends in a sensitive way*on the 
electron configuration and, in transition metal alloys in 
particular, the d shell electrons are thought to conten 
Sronificantly. | 

Very little work has been done on the elastic constants 
of binary alloys where both constituents are transition 
metals. In addition to making an isotropic alloy some 
knowledge about the interaction between atoms in such an 
alloy would be expected to be gained from the analysis of a 
series of compositions. Therefore, seven Salleerstelil Aine 
compositions of the vanadium-iron system were produced with 
Miron content between 3.25 and 15.32 at.%. The longitudinal 
and transverse elastic constants were determined, by a pulse 
echo method, for each composition. Then, using an empirical 
relation for the shear constant Cisi/2(C -C, i), the single 
crystal elastic constants were calculated and discussed in 
light of the variations in electron to atom ratio and the 
Mewbanad Contributions. 

The classical theory of elasticity is discussed in 


Appendix I-and the Fuch's theory of determination of the 





elastic constants from fundamental considerations is discussed 


in Appendix II. 





II. EXPERIMENTAL DETAILS 


A. SAMPLE PREPARATION 

The preparation of the seven alloys used in this study 
was done at the Naval Research Laboratory, Washa nese Deer 
The history of each sample measured is shown in Table I. 
Commercial analysis of the composition and impurities of 
each sample is shown in Table II. 

Two methods of producing the alloy were used. The 
constituents were melted by an electron beam in vacuum or 
arc-melting under an ultra-pure helium atmosphere. Electron 
beam melting was discontinued after the second sample when 
it was determined that the refining effect of this method 
reduced the iron content of the alloy by unpredictable 
amounts. Ingots from both methods were hot forged and 
Swaged into rods. A typical rod was 3/8 to 1/2 inch in 
diameter and 12 to 15 inches long. 

tmese rods awere used ln attempts to guow alloy siege 
crystals. The attempted growth methods were liquid float 
zoning and strain annealing. All of the samples exhibited 
grain growth: zone melting produced grains of .25 to .75 
mm diameter, while strain anneal methods produced grains 
of 1.25 to 3.75 mm diameter. Time and the availability of 
material precluded a more thorough exploration of the pos- 
Sibilities of producing single crystals and thus poly- 


crystalline rods of each composition were selected for 


measurement. 





Fourteen specimens were made from the seven alloys. 
Seven were 2 cm long and .635 cm in diameter. These HOLE 
used for measuring the longitudinal velocity. The remain- 
der, cut adjacent to the first, were 1 cm long and were 
used for measuring the transverse velocity. To provide 
the best possible diameter to wavelength ratio in the 
transverse samples the largest available diameter was pre- 
served in each case. Sample diameters varied from .635 cm 
to 1.10 cm. The ends of all of the cylindrical specimens 
were machined flat and parallel and the surface polished. 
The surfaces were found to be parallel to an average of 


1.5 milliradians. 


B. MEASUREMENT OF THE SOUND VELOCITY 

The velocity of sound in the seven alloys was measured 
by the pulse echo technique. A sound pulse was transmitted 
into the material by a piezo-electric quartz transducer 
driven by a pulse generator operating at the transducer 
iinmcdanental frequency or an odd harmonic. ©The initial 
pulse and each succeeding echo can be received by a second 
transducer. This second transducer is placed at the oppo= 
Site end of the sample. For the transverse measurement the 
two transducers must be aligned so that their particle 
motions are parallel. The two transducer method eliminated 
saturation of the receiver amplifier by the driving pulse. 
The arrival time of each echo is usually determined by 


using an oscilloscope trace with a delay time multiplier. 





The transducers used were 1OMHz X cut quartz for the 
longitudinal mode measurements and AC cuts of the same 
frequency for the transverse mode. The transducers were 
bonded to the samples with Salol. A spring loaded jig was 
used so that both transducers could be bonded simultaneously. 
This improved the uniformity of the bonds and eliminated 
the problem of melting the first bond while the sample was 
being heated for the second. During the measurements the 
sample was suspended at its center by a three point circular 
holder ee concentric with two shielded contacts spring 
loaded to the transducers. Series coils were employed to 
provide a low impedance match of the capacitive transducers 
to the pulse generator and receiver. All measurements were 
performed at ambient room temperature (294 K). 

The alloys in this experiment were, in most cases, highly 
attenuating. The transverse mode was generally more at- 
tenuated than the longitudinal mode and often only two or 
three echos could be discerned above the system noise. 
Measurements of the total attenuation were made in the 
longitudinal case. The total attenuation includes all loses, 
which are loses at the bonds and transducers as well as in 
the sample themselves. These loses varied from a low of 
3db per cm for the best sample to approximately 5 db per cm 
for the worst sample. For this reason the unrectified pulse 


waveform, after amplification, was further processed by a 


La) 





Princeton Applied Research (PAR) Model 160 Boxcar 
INnveqracea. 

The boxcar integrator has three principal advantages for 
the pulse echo experiment. It provides a large improvement 
in signal to noise ratio. It faithfully reproduces the small 
amplitude portions of the waveform so that the arrival time 
of the pulse can be accurately determined. In conjunction 
with a strip chart recorder it provides a permanent working 
data record. The boxcar integrator operates by scanning an 
aperture in time along the waveform and averaging the value 
of the signal at each point over many repititions of the 
waveform. The result is that the random noise fluctuations 
on the signal will average to zero, and if enough repetitions 
of the signal are integrated, the output will be the true 
value. Typically 50 usec of the echo train was scanned 
using a 10 nsec aperture. Such an aperture can resolve 
frequencies up to about 30 MHz. The scanning time to make 
each pass along the train was 50 minutes. A pulse repitition 
frequency of 10 KHz was used so that each point on the wave- 
form was sampled 6xl10° times. The output to the strip re- 
corder is then a picture of the echo train produced at one 
microsecond per minute. By diverting the boxcar input to 
the time mark generator that was triggering the experiment 
a trace of convenient time marks could be placed in the 
"dead space" between echos. These time marks could not be 
used to measure the arrival time of echos, however, since 


any variation in the linearity of the travel of the aperture 


del 





along the signal (ie along the timebase) or any variation 
of the speed of the strip chart recorder would have the 
effect of expanding or contracting the trace. Instead the 
voltage corresponding to the position of the aperture along 
the timebase was used to drive a second pen which ran con- 
currently with the echo trace. (Figure 1) This second 
trace was used to determine the echo arrival times. The 
accuracy of the time measurement 1S improved by this re- 
finement since if the nonlinearities above occur they will 
be automatically compensated. The time marks on the echo 
trace were used to check the calibration of the time trace. 
Electronic instability over the relatively long times re- 
quired to make a trace did not present a problem. Varia- 
tions in the trigger rate of the time mark generator or in 
the trigger response of the pulse generator and boxcar 
integrator that would have resulted in noise in the output 
beace were not observed to oeeur. From these traces the 
arrival time of each echo was measured and plotted against 
the elapsed path length. The slope of this plot is the 
sound velocity. 

In a homogeneous polycrystalline material the value of 
the transverse velocity should not depend on the orientation 
of the transducer with respect to the sample. Such is not 
the case in single crystal measurements where alignment of 
the particle motion with certain crystallographic directions 
1s required to determine the elastic constants. The rather 


large grain size in these samples could have the effect that 
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the orientations of the grains might not beysutfiverencly 
random to result in the true average velocity being measured. 
This effect of a lack of randomness is reduced if the aniso- 
tropy of the individual grains is low, and vice versa. To 
compensate for any lack of random grain orientation velocity 
measurements were made on each sample with the induced 
particle motion oriented at 0, 60 and 120 degrees from an 
arbitrary mark. The velocities thus measured and several 
values at unknown orientations were combined by taking the 
arithmetic mean as the velocity for the composition. The 
choice of a difference in angle of II/3 for the various 
orientations was chosen as convenient and because a cubic 
material has II/, symmetry. The variation of sound velocity 
with orientation for each sample is indicated by the bars 

em the data in Figure 4. These bars are drawn between the 
maximum and minimum velocities observed. The average varia- 
tion in velocity between orientations and samples was 2.26%. 
The repeatability on a given sample for a given orientation 
is within the accuracy of the experiment which is discussed 
below. 

To determine the elastic constants both the sound 
velocity and the density of the sample must be known. 
Pearson [Ref. 2] lists data by Hanneman and Mariano which 
includes the value of the lattice parameter of the vanadiun- 
eon system for 0, 9.5 and 20.5 “ates iron. “Ltva-sneocn 
curve is drawn through these points ene lattice parameter of 


the seven alloy compositions used can be found. The density 


a3 





determined by this method is compared in Figure 2 with the 
density determined by weighing and measuring the samples. 
Figure 2 shows density as a function of the atomic percent- 
age of iron and indicates the agreement between the two 
methods is within 1% for all compositions. The experimental 
points can also be fitted to a straight line with one point 
outside of 1% and an average deviation of 0.3%. Such a 
straight line fit would increase the density determination 

of the dilute composition approximately 0.5% and decrease the 
concentrated composition a like amount. It is felt that the 
precision of the x-ray determination of the lattice parameter 
exceeded that of the experimental check on the calculated 
density so that the density of each composition has been 
calculated from the lattice parameter read from a smooth 


curve through the published values. 
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Tit. EXPERIMENTAL RESUS 


The elastic constants of a cubic polycrystalline material 
are the shear constant G, and the longitudinal constant, 
(B + _ G), where B is the bulk modulus. Appendix I shows 
that the longitudinal and transverse sound velocities are 


given by 





In order to compare the quantities measured for the 
alloys with those of pure vanadium the Hill method [Ref. 3] 
was used to average over all lattice orientations the single 
crystal elastic constants of vanadium as measured by Alers 
[Ref. 4]. This average is required to obtain the poly- 
crystalline value of G and is also discussed below. The 
bulk modulus in cubic materials is the same for the single 
crystal and polycrystal case. After a small correction 
for the difference between the density used by Alers and 
that calculated from the lattice parameter B and G for pure 
vanadium are known. For vanadium and each composition, the 
atomic percentage of iron, density, longitudinal constant, 
shear constant and bulk modulus are given in Table pl a 
Graphs of the measured elastic constants which also contain 
the calculated values for vanadium are shown in Figure 3 


to Figure’ 5. 


i> 





The longitudinal modulus (Figure 3) increases with in- 
creasing iron concentration. The experimental data is fitted 
best by a least squares parabola. With the exception of the 
3.25 at.% composition, which is approximately 2% low, all 
of the other samples and the vanadium value fall within 13% 
of this fit. Above approximately 5 at.% iron the measured 
values increase essentially linearly. However, a linear 
least squares fit to the data had a larger least square 
error and fell approximately 3% low for the value of pure 
vanadium. 

The transverse modulus (Figure 4) is essentially con- 
stant over the entire range of compositions investigated. 
These is considerably more scatter in the data for the 
transverse case than in the longitudinal case. The best 
fit straight line increases approximately 1.5% over the 
pure vanadium value for the 15.32 at.% composition. All 
of the samples fall within 5% of this line with the ex- 
ception of the 7.40 at.% composition which is approximately 
12% low. The reproducibility of this low velocity for 
sample 3 was demonstrated on a third specimen cut adjacent 
to the first two. Because of the small size of the transverse 
constant with respect to the longitudinal constant the 
percentage errors are larger although the absolute errors 
are approximately the same. 

The bulk modulus (Figure 5) is readily obtained from 
the definition of the longitudinal Be rae eh 5 G). 


Since it is a combination of both the longitudinal and 


mG 





transverse constants it reflects both the rise with in- 
creasing iron concentration of the longitudinal constant 
and the scatter present in the transverse data. 

Based on the estimate of the errors in the measured 
time, length and density the accuracy of a single measure- 
ment of ov is estimated to be approximately .75%. The 
largest portion of this error is in the time decree 
Differences in individual echo arrival times for different 
runs on the same sample that resulted in the same velocity 


were found to be approximately .75%. 
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IV. DISCUSSION 


The two elastic constants measured in a polycrystalline 
experiment are the shear constant, G, and the longitudinal 
constant (B + _ G). For this experiment therefore the 
experimental results are the variation of the shear modulus, 
ee eana the bulk modulus, "8, with increasing atomic per— 
centages of iron. 

Although the interpretation of these results alone is 
interesting they do not provide as much information, that 
would be useful in interpreting the fundamental contributions 
of atomic forces to the elastic constants, as the single 
crystal elastic constants would. (See Appendix II for 
the theory of the fundamental contributions to elastic 
constants) Moreover, for cubic materials there are three 
non-zero independent elastic constants, (See Appendix I) 
and it iS impossible to unfold these three quantities 
from the two measured quantities. 

An independent source of the third constant was avail- 
able in the work of Fisher and Dever on the stability of 
bcc transition metals [Ref. 5]. They present an empirical 
relation for the C' elastic shear constant as a function 
of the electron to atom ratio of transition metal alloys. 
Here C' is one of a set of 3 constants which are linear 
combinations of the basic cubic constants, Ci: Ci2 and 


Cyuy (See equation Al.8,f.f.). 
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B= 
3 
erie Oe (1) 
C'= ——____ 
2 
C = Cyy 


Here B is the bulk modulus and C and C‘' are shear constants. 
In terms of this set of constants the measured bulk modulus 
for the polycrystalline material is directly comparable with 
B=1/3(C11+2Ci2), however, the shear modulus, G, is a function 
Se poth C and Cc’. The, fact that Gis a known function er 

C and C' is a result of the methods of determining poly- 
crystalline aggregate moduli, as averages of the single 
@evstal moduli. It is thus instructive to digress bricily 

to discuss these averages. 

In 1952 Hill [Ref. 3] showed that two classical approxi- 
mation to polycrystalline averages by Voigt in 1928 [Ref. 6] 
and Reuss in 1929 [Ref. 7] were upper and lower bounds on 
the problem and suggested the arithmetic mean as an average 
that would be closer to measured values. Voigt averaged 
over all lattice orientations using the assumption that 
the strain is uniform throughout a grain. Reuss did the 
averaging assuming that the stress is uniform throughout 
a grain. Following Hill, the strain energy function of 
an isotropic material in terms of the principal stresses 
T. and strains Be and the unknown bulk and shear modulus, 


B and G are 


Wg 





2 
w = B(e)+e,+e3) 7+ -G Janie -es). (Dy 
3 i¥j=) an 
3 
rar Hn CWO oa 3 
W= —(T,+T4+T Se Toole 
B Pod) 5G teen ees (3) 


Even though the Voigt and Reuss approximation are limits 
they lead to simplified relations between the polycrystal- 


line constants and the single crystal constants. 


9B = (Gree oC 345) 24 C7) o+Co 3+C3 7) (4) 
eGo =) (CritCoetC 33) =eeGie+eo44C31) =F3 (CygtesteCcem (>) 
- = (S11tS224S33) + 2(Si2tS23+Ss1) (6) 

r 
aa SCS pitosstS33) SAiS ot So3tSai to (Sua toceaoe (7) 

06 


These relations apply to all crystal classes and the sub- 
scripts v and r deonte Voigt and Reuss respectively. The 
244 are the elastic compliances (see Appendix I). 

5 OF Sis in (4)through (7) in terms of 
mre cubic constants C)1, Ci2 or Cy, and using the cubic 


By writing the Cc. 


miverse relations 


il 


eee (Cee 20>) 
(Siit+2S512) 


(S11-Si2) 


Cyy= (Cig Ci) 


abs 
Sun’ 


one obtains, 


1 
B = B = = ({C71t2Cj19). (8) 


Z0 





which is identical to the single crystal experssion in 


equation (1), and 


: Fenn (GS 
G= Z(Cri-Ci2t3Cus) G = SCyy(Cii-~Ci2) . (9) 


M AC yo (CnC) 


Now converting to C and C' using equation (1) 


2 3 ' 
aoe - CO g¢ = =o __.. ALO) 


4c + 6C' 


The Hill average of these G = 1/2 (G_ ats G..) is 


1 petere > + 5.2cc' 
4c + 6C! 


G= Ci) 


Equation (ll) is the desired expression for the shear 
constant in terms of C and C'. 
Returning now to the empirical relation of Fisher and 


Dever for the shear constant C', 
c' = 0.01391 x 10?? (e*/a) 3°? *ayne/cm?. (12) 


This power law relation for C' was determined by taking 
the best fit of values of C' for three Ti-Cr alloys and the 
published values of C' of vanadium and chromium. In this 
fit the electron to atom ratio was a parameter. The details 
will be considered below. The central result was that by 
assigning effective electron to atom ratios, e*/a, of 2, 
3, and 4 to titanium, vanadium and chromium respectively 
they could obtain the relation (12). 

Two important cuestions can be raised. First why does 


C" seem to depenc more on an electron to atom ratio, than 


eal 





on the electrostatic and ionic contributions to the shear 
constants outlined in Appendix II. Second, why the par- 
ticular numbers 2,3,4 instead of the total valency or some 
other number. Quantitative answers are not forthcoming, 
however qualitatively the reasoning is as follows. If in 
the region of (e*/a) from 2.15 to 4.0,C' is governed by the 
relation (12) (the e*/a range of the present work is from 
3.0 to 3.46) then a relation of the same form and similar 
magnitude should exist for the 4d and 5d transition metals, 
Nb, Ta, Mo and W. By plotting the published values of C 
c' and C,versus e*/a for -V,; Nb Ta, Cr, Mo, and wand looming 
at the percentage variations between the 3d, 4d and 5d 
elements they conclude that although C' does not depend 
entirely on e*/a it has much more of this type dependence 
mmam dO C;}1 or C. mheyespeculate thatewdrelecerenme.— 
change interaction with nearest neighbors could overcome the 
negative contribution of ionic overlap discussed in Appendix 
er. 

For the choice of e*/a of 2,3 and 4 for titanium, vanadium 
and chromium respectively Fisher and Dever refer to a 
recent APW band calculation of the transition elements by 
Snow and Waber [Ref. 8]. In this calculation, Snow and 
fewer CcOlpure the electron occupancy of the 5, p and a 
bands as well as the number of electrons in the external 
region between the atomic spheres that are characteristic 
of the APW method. These SUT ci SE computed for two 


Mors 


Rene. . n 
initial configurations of electrons, s‘d and s*d where 
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n varies from 2 for titanium to 9 for copper. Snow and 
Waber favor the gigh*1 configuration as the most probable 
and, for this configuration, show d shell occupancies of 
2.11, 3.07 and 4.01 for titanium, vanadium and chromium 
respectively. These are very close to n which is the 

e*/a chosen by Fisher and Dever for the same elements. Thus 
if the d shell interactions are considered the primary 
contribution to C' this seems to be reasonably good evidence 
for the chosen values of e*/a. 

For the vanadium-iron alloy of this work the e*/a of the 
seven compositions could be computed using the value of 
3.0 for vanadium and by referring to Snow and Waber's work 
choosing the ratio 6.0 for iron. They give ad shell oc- 
cupancy of 6.18 for iron starting from the stat! configura- 
tion. Then using equation (12) C' could be computed for 
each composition. These values of C' and the measured 
values of G can then be used to solve equation (11), for 
the shear constant C of each composition. These results 
are listed in Table III along with the best fit values of 
the derived quantities. 

In Figure 6 is shown the results of all the known shear 
eonstant data for the bcc region of the first period tran- 
Sition elements and their alloys. The values of C and C' 
shown for vanadium-iron region are speculative in that they 
are derived from an empirical data fit and a polycrystalline 
experiment. However, the derived values of C in the reqion 


Sy a = 390 Go 3-46 are continuous with the results creme 
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Fisher and Dever alloy experiment which are in the region 
e*/a = 2.15 to 2.6 of the figure. The only other known 
experiment in this region is the work on molybdenum-rhenium 
by Davidson and Brotzen [Ref. 9] which is shown in the region 
ex¥/a = 4.0 to 4.26 of the figure. Fisher and Dever had no 
explanation for how the values of C for vanadium and chromium 
joined up. These results of the present work merely state 
that C decreases before increasing to the chromium value. 

The scatter in the C values reflects the scatter in the 
experimental results. The points for V-Fe alloys on the C' 
line in Figure 6 are calculated from the empirical relation 
(12) and it can be seen that C' apparently has a maximum 

at e*/a = 4.0. Fisher and Dever point out that this maximum 
corresponds to the maximum d band width in Snow and Wabers 
Galculation. 

Variation in the anisotropy ratio with e*/a is readily 
feparent in Figure 6. At the left hand sade the atic ds 
greater than unity and decreases to that value for e*/a 
of 2.8. It continues to decrease through the V-Fe region. 
In the region above e*/a = 4 the anisotropy is again in- 
creasing for the Mo-Re alloys. The anisotropy returns to 
unity for e*/a = 4.26. Thus the concept of an isotropic 
V-Fe alloy is probably incorrect. A more likely choice 
would be vanadium with a few atomic percent titanium to 
make e*/fa = 2.8. Ti-V forms a bcc solid solution in this 


region. 
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It is interesting to note in Figure 6 that C extrapolates 
to zero at an e*/a of approximately 4.2 which also heralds 
the onset of the sigma phase of the vanadium-iron system 
fRef. 10]. That C would continue to zero and create an in- 
stability seems unlikely and it is thus assumed that C 
must increase eventually toward the Cr value. In terms 
of the energy of an electron gas discussed in Appendix II, 
the addition of electrons in the region above e*/a=3.5 would 
push the Fermi surface closer to the Brillouin zone boundary 
where shear strain on the direct lattice, which is, of 
course, reflected in the reciprocal lattice could increase 
the average electron kinetic energy and contribute a positive 
component to C. Another way of visualizing the electron 
gas effect on C is in terms of the density of states function. 
Snow and Waber.show that the density of states at the Fermi 
level varies with group number in such a way that a minimum 
occurs for e*/a just less than 4 for the gignht} COnmietiiation . 
This is the region where increasing e*/a causes the Fermi 
level to traverse between the lower energy bonding states 
of the d band and the upper energy anti-bonding states. The 
correlation between these effects and the requirement that 
the C value reach a minimum and turn up toward chromium is 
unknown. It 1s only pointed out that they both occur in 
the same region of e*/a. This concludes the discussion of 
the shear constants. Before going on to the bulk modulus 
a short discussion of how the various Pease. e*/a, were 


determined is included. 
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In making the choice of e*/a for the elements, Fisher 
and Dever plotted C' for their three alloy compositions, 
vanadium and chromium versus l, 2, 3 and 4 excess electrons 
in chromium over those in titanium. Only the value 2 gave 
a smooth curve through the five points. They then plotted 
mo C' versws ln e*/a for and e*/a of Cr from 1 to 8. "The 
only e*/a to come within the experimental error of C' for 
their alloy and V was e*/a = 4 for chromium. This fit was 


oye 5 


oo 


low for the published C' of Cr. and had the parameters 
given in (12). 

Since they had established the e*/a for V the remaining 
choice in this work was the value for Fe. Values of C were 
calculated using e*/a = 6 and 7 for iron. An e*/a of 7 
has the effect of increasing the maximum e*/a for V-Fe to 
3.61 and slightly lowers the C values. However, since G 
is very insensitive to change in e*/a a choice between 6 
and 7 could not be made on this basis and it was decided to 


Pe@consistent and use the s!a™™ 


configuration with an e*/a 
of 6 for iron. If possible errors in the determination of 
G are considered it is found that even if G were to increase 
with alloying at the same rate as the longitudinal constant 
the physical result that C decreases with alloying for e*/a 
above 3.0 would only be changed in degree. 

Published values of the bulk modulus in the transition 
metals [Ref. 11] show that it is peaked in each long period 


in the vicinity of Group VIIIB. This dependence is similar 


to that of the cohesive energy and is independent of an 
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interpretation of any model of the contributions to the 
bulk modulus. Converting to electron concentration, as has 
been discussed above, these data and other published results 
for alloys support the idea that the bulk modulus increases 
with increasing electron concentration up to a maximum at 
ex*/a of 4 or 5. 

One may think of the buik modulus as arising from the 
"free" electron gas and from the ion-ion interactions. 
Snow and Waber using the giant? configuration estimate that 
the s like portion together with the charge external to the 
atomic sphere is between 1.0 and 1.5 electrons. [In the 
alloys investigated here the lattice shrinks with the 
addition of solute. This means that if only these electrons 
are considered as contributing to the bulk modulus it will 
increase with decreasing lattice parameter. A simodle 
estimate shows that this contribution is about 30% of the 
total for pure vanadium and roughly the same fraction for 
the alloys. ‘Thus one may infer that both contributions 
increase with alloying. Then for the ionic contribution 
one may conclude that alloying iron with vanadium increases 
the bulk modulus either because the iron-vanadium ion inter- 
actions are larger or because the iron-iron interactions 
are larger. From the shape of the bulk modulus curve the 


latter interpretation seems more probable. 
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VV, CONCLUSIONS 


Addition of iron to vanadium increases the bulk modulus 
but not the shear modulus, although the latter data are 
scattered. The increase in the bulk modulus can be partially 
attributed to an increase in electron gas concentration and 
partially attributed to an increase in ion-ion interactions. 
The relative contributions of the s and "free" electrons 
maceugh the electron gas and 1on-ion exchange interaction 
is approximately constant for the range of alloys studied. 

Since the shear modulus remains constant this means that 
the shear constants C and C' do not change with alloying or 
more probably that they each change at different rates. 
Using the empirical relation of Fisher and Dever for C' 
the variation of both shear constants with electron concent- 
ration is consistent with the previously measured values 
for vanadium and the titanium-chromium alloys. 

The anisotropy ratio is a smoothly varying function over 
the range of effective electron concentrations from 2 to 
4.5. Thus this work supports the idea of an sight} con- 
figuration for the transition metals and an almost complete 
dependence of C' on the electron to atom ratio. The concept 
of an isotropic alloy of the transition elements has been 
demonstrated in the Mo-Re case. The existence of a second 
such alloy between two first long period elements, V-Ti, 


is indicated. 
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Future work in this area could include investigation of 
the elastic constants of a series of various alloys with 
e*/a from 2.15 to 4.30. Such alloys as V-Ti, V-Cr, V-Fe 
and Cr-Fe might be employed to investigate this entire 


region of bcc stability in the transition elements. 
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No. Production 


technique 
ae arc melting 
2 electron 


beam melting 


3 arc melting 
4 arc melting 
5 Sleceron 


beam melting 


6 arc melting 
i arc melting 
Notes 


TABLE I 


Vanadium 
source 


Argonne 
National 
Laboratory 


Bureau of Mines 
electrolytically 
refined 


Bureau of Mines 
vacuum arc 
melted 


Bureau of Mines 
electrolytically 
refined 


Naval Research 
Laboratory 


Argonne 
National 
Laboratory 


Bureau of Mines 
electrolytically 
refined 


Crystal growth and 
anneal technique 


annealed 11/2 hours 
at 1000C longitudinal 
and transverse velo- 
city were measured 
annealed 4 hours at 
1600C 


Float zone melting 


strain anneal 
annealed 21/2 hours 
at 1600C 

swedged to strain 
annealed 11/2 hours 
at 1600C 


passed slowly through 
an ianduct#on coud 
after straining 

this took 4 hours 

at 1600C 


float zone melting 


Same as sample l 


float zone melting 


1. Iron for all samples was produced at NRL [{Ref. 12] 


2. Vanadium produced by the Bureau of Mines is reported in 


(REEwelS angela). 


3. All anneals were done in vacuum 


4. Temperatures during anneals were determined by optical 


Dy POMeeLy . 
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V 0.00 6 
i B22) ae 
2 Ds 35 EO 
5 7.24 26 
3 7.40 6 
4 10.96 6 
bee 13.94 6 
fol. 352 6 
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5) 
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9 in gm/cm? 
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OMe), 
18309 
.2438 
2209 
3 01/7 
16 
- 4609 
~4917 


No NM N NO NM NY NY NH 


SQ © 02 0 OQ Oo Oo 
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~474 
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~452 
-490 
~415 
~479 
-483 
-469 


TABLE III 


e*/a Cc C B 
Bevo 0.546 0.431 Pe So4no 
BevlOmeONG09 40.410 1.5 
S16 90.1649 0.353. 1.6Z2eun0 
3.22 02691. 0 cee ele scs eo 
3222 OR69l SOe2 oie scree 0 
3233 027773 ~0.344 lease 0 
3242 02845 05526 ) tae Oo 
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APPENDIX £ 


The elastic theory of solids is based on the Hooke's 
Law assumption that for sufficiently small strains, SE 
stress, T. is proportional to strain, where the elastic 


constants, a9 are the constants of proportionality. 
6 
cs = os. tere. Alaa 


Deformation of the medium in an adiabatic manner results 
in a stored elastic energy density, expressed by a strain- 


energy function, w. 
6 
EWG aes ces Al.2 


In this notation the relation between 1 or j and the co- 


magnate direction is 


Here 1, 2 and 3 represent the normal strains which are 
those resulting from a stress on the face of a cube normal 
to the x axis and in the x direction etc. The numbers 4, 
5 and 6 denote the shear strains which are those resulting 
from a stress on the face of a cube normal to the y axis 
and in the z direction etc. The shear strains in terms 

of infinitesimal displacements are set in equation 
Al.15. In both cases the summations on i and j allow for 


the general case of an anisotropic medium. 


Ss, 





In terms of the energy density and/or the stress the 
constants of proportionality, ory called the elastic 


constants, stiffness constants, or moduli of elasticity 


are 
aT. ; 
= AE a Ob eW. 
ij de.  2e. oo (Al . 3) 
J a J 


It is sometimes convenient to express relationships in 


terms of the reciprocal of the set of equations (Al.1). 
6 
e. = 2 §.. T. (Al .4) 


The coefficients of this set of equations, yy are called 
the moduli of compliance. Treated as matrices the product 
matrix of the one! and ic.) array is the unit matrix, 
[S53] oa = I. 

The strain-energy function, w, is the work done per unit 
volume in deforming the crystal and thus has the dimensions 
of erg-cm*®. The strains, e, are dimensionless, thus the 
elastic constants also have the dimensions of an energy 
density usually expressed as dyne-cm °. Inebhl Sewer cede 
has been convenient to show the C,. in units of 10+? dyne- 
cup. 

The number of independent elastic constants is not 36 


as equation (Al.2) might imply. From equation (Al.3) and the 


fact that é6w is a perfect differential we see that 


“4-7 0GT = semdeaa semoen > xe San (Al.5) 
J 1 J 7 1 
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This symmetry See reduces the maximum number of indepen- 


Cia 
ji 

dent elastic constants to 21. The elastic constants are 

usually shown ina 6 by 6 array which for the most general 


@ase 1s 


iC 1 2 gis Cupane 156 < 1c 
Gite C72 €o 3 ©5777 e>5 Cc 
Cryswe> 3. C337 es wees Cac 
(AlL.6) 
Cire 2 hate Cyn Crs. Cig 
Cys eo 51°Cs 5. Cys Cs6 Caz 


Cié ¢c26 C36 Cue C56 Cee 


7 


The number of non-zero independent elastic constants is 
further reduced, except for the triclinic lattice, by the 
effect of lattice point symmetry operations. This effect 
can be demonstrated by writing down the strain-energy 
function in an original coordinate system and one obtained 
from the original by a lattice symmetry operation. If the 
two strain-energy functions are equated and coefficients 
of like strains compared the results are that certain 
constants become zero and others equal. For the particular 
case of cubic symmetry only 3 constants remain Ci1, Ci2, 


and Cy,. Thus the array (Al.6) becomes 


Al 








Ci2 Cet Ci2 0 0 0 

Ci2 Ci2 Cii O 0 0 
(Ae) 

0 0 O Cus 0 0 


Using this array we can write out equation (A1.2) for the 
strain-energy function in a cubic crystal. 


1 1 


Meee. [le te-+e-l+Cy>olemeate G@ te e {4-€ 9 fers 
2 lees coon ame eee 3. 3 


& S) 6 


(A1l.8) 


A direct application of equation (Al.3) defines the cubic 
elastic constants ie C = a etc. It is often more 
convenient when dealing with theories that derive the elastic 
constants from fundamental concepts to define certain 


combinations of the constants C1, C2 and Cyn. AN equiva- 


lent set of these constants are 


Ge = 1/72. (Ci q-Er) 
C = Cys 
B = 1/3(Cii+ Z Cio) (Al .9) 


In this set the constants C and C' are shear constants 
chosen because they can be obtained from the strain energy 
function by strains that preserve the volume of the unit 


Geil to first order in the e.. The constant B, is the 
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bulk modulus and is the elastic constant corresponding 

to hydrostatic compression. These three cases are ob- 

tained from the cubic strain-energy function (Al.8) by choice 
of the appropriate strains. 

1. C' = 1/2(Ci11-C}2). %In this case the cube is compressed 
along one axis and extended along a second so that two of 

the strains, Say, ey and e, are equal in magnitude but 
opposite in sign. All other strains are zero. Then 


equation (Al1.8) becomes 


w= (Ci1-Ci2)e17 from which 
C'= L a*w | ESI) 
l de, cee . 
2. C= Cyy. This corresponds to shear ina plane parallel 


to a cubeface so that one strain, say e, is finite and all 


others are zero. Then equation (Al.8) becomes 


A) eC eyrerre from which 
_ o¢w 
Se dey | aye ee 


See B = 1/3(Ci1 + 2 Ci2) For the case of hydrostatic 
compression the three normal strains are equal and negative 
and the three shear strains are zero. Then equations (Al1.8) 


becomes 


w= = (ORS ees = from which 
Bee 37w 
9 J0@}] 2; = Oo (A1.12) 
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A cubic crystal is said to be isotropic if the two shear 
constants in equation (Al1.9) are equal. The ratio of these 
constants, A = C/C' is called the anisotropy and for an 
isotropic crystal A=1l. Thus, the effect is to reduce the 
number of non-zero independent elastic constants to two. 
These are generally chosen to be Cyy = uw and C)2 = A. 

The constants A and u are called the Lame' constants. 

The elastic constant C;, then maybe written in terms of 

the Lame' constants and the elastic constant array (Al.7) 


contains only two independent constants. 


(A+2u) r r 0 0 0 
r (At2u) A 0 0 0 
r A (At2u) O 0 0 
(Fale 13) 
0 0 0 u 0 0 
0 0 0 0 u 0 
0 0 0 0 0 u 


meaaeequation (Al.1) the 1soetropie stress-strain relacronsaimay 


now be written out. 


T, = (A+2u)e,+ Aleote3) = AAt2ue, 


T2 = (A+2u)e2+ Al(eite3) = AA+2ue2 
T3 = (A+2uU)e3+ Ale2te;,) = AAt2uUe3 
(Ade i4) 
Ty = yes 
mee Ves 
Te = wes 
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Here A = e,+e2+e3 18 called the dilatation. Following Mason 
[Ref. 15] the expressions for the transmission of sound waves 
in isotropic solids may now be readily developed. 

In terms of infinitesimal particle displacements u, v, 


and w in the x, y and z directions respectively the strains 


are: 
dU OV OW 
e)= <Guee So — —— e 3= es 
OX oy OZ (A1.15) 
ee pee CMe ee (ov _ ou 
i dZ«OOY : dZ «OX 5 ax oy 


The equation of motion of an infinitesimal volume element 


of density in the x direction is 


2 - 
fe) “4 dx dy dz = ie (Alm ) 
meeesimitarly for y and z. “me force F, is that arising 


from a stress -T,(x) on the face at x of a cube of volume 
AxAyAz and a stress T; (x+tAx) = T; (x) + Ax on the face 
pee xXtAx). Other forces in the x direction arise from 


the shear stresses T, and 1 so that in equation (Al .16) 


OT ¢ 
OX dy 


+ So8)dx dy dz UML) 
and similarly for Mo and Fo 
It one substitutes for the e. in equation (Al.14) 


equation (Al.15) and take the derivatives indicated by 


equation (Al.17) the equations of motion (Al.16) become 
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Sail oA m 
p jt 2 oa a2 a) aX Stoemt Yamal 
ay o, Tage 2 
0 a2 = (A+2 u) ae uV “Vv (Ade Ss) 
Saw dA 2 
Oayee = (A421) 7S uv “w 
2. 9? NG AVE 
petro |6 = ee ay? ane? 


If the three equations (Al1.18) are differentiated with respect 
moe, Y and z respectively and the resulting equations 


summed, a wave equation for the dilatation is obtained 
soe Ne A, 2 
aes (A+ 2) VA (A1.19) 


Por; a plane wave with partic¥e motion u only equation (Al. 19) 


becomes 
ae ule 97u : 
O ate (A+2,u) apres (A1l.20) 


which shows the velocity of a longitudinal wave traveling 


tmaene x direction to be 


_ A+24 
Ve = pam (Al.21) 


To obtain a transverse wave equation the first and third 
of equation (Al1.18) are differentiated with respect to z and 


x respectively and the difference is formed. 


9° JW owes => cou Weow 
Payee oe gale Uae oe (Al. 22) 


Then assuming a transverse plane wave traveling in the z 


direction with particle motion in the wu direction, 
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= gn WSEAS) equation (22) becomes 


Co d7u 


from which the transverse wave velocity is 


v = {ft | (Al.24) 
S , 


Finally, in the discussion of polycrystalline materials 
whose isotropy results from the random orientation of 
generally anisotropic crystallites the shear constant has 
been called G. Setting G equal to uw and using the definition 


of the bulk modulus along with the Lame' constants 
(A+2u) = (B + A G) (A1.25) 


Seostituting into (Al.21) amd (Al.24) the polycrystalline cubic 


material velocity equations become 


Vy, = 





(Al.26) 


47 





APPENDIX II 


The treatment of elasticity in Appendix I is based on a 
classical Hooke's Law force and its attendent strain-energy 
function. In calculating the elastic constants from funda- 
mental consideration one must generate a function U which 
is the total energy per atom. The equivilance of this 
function and the strain-energy function, w, iS assumed so 
that in terms of the energy, U, the definition of the 


elastic constants is 
(A2.1) 


The total energy is a sum of contributions arising from 
the various forces between atoms in solids. 

A successful theory of U which is a function of the 
interatomic spacing will give first, the cohesive energy 
of the crystal and secondly upon differention, the equilibrium 
lattice constant. Thirdly, another differentiation of U 
as in (A2.1) will give the elastic constants. This latter 
mesULtrement 1S @ rigorous test of ‘the validity of any 
fundamental theory of the total energy. 

Present day theories of the total energy U, which have 
been applied to the determination of the elastic shear 
constants C and C', originated with Fuchs. Fuchs applied 
the methods of Wigner and Seitz for determining the cohesive 


energy of metalic sodium. He extended the treatment to 
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copper, and later calculated the elastic constants of 

eepper (Ref. 1)< 

| Fuch's calculation was applied to transition elements 

by Zener [Ref. 16] and Isenberg [Ref. 17]. They concluded 
that the observed stability of the bcc transition elements 
resulted largely from a repulsive interaction between second 
nearest neighbors arranged in a anti-ferromagnetic array. 
Neighbors and Smith use Fuch's calculation to discuss dilute 
alloying of other elements with copper [Ref. 18]. 

It is instructive to consider the calculation of Wigner 
and Seitz for the cohesive energy of metalic sodium since 
the contributions to the total energy for this case can be 
extended to other monovalent metals. 

The electronic configuration of sodium is a single s 
electron outside of a closed neon core. In the solid 
state this electron is assumed to be free to move through 
the crystal lattice and Wigner and Seitz calculate the energy 
of the lowest state, Uss of the free electron. The cal- 
culation is facilitated by imagining the lattice to be 
divided into polyhedra constructed by erecting planes which 
are perpendicular bisectors of the lines joining an atom 
to its neighbors. Such polyhedra are adjacent and fill 
all space occupied by the lattice. The motion of the 
electron through the lattice then becomes a motion between 
adjacent identical polyhedra. Such a motion will be period- 
ic with respect to the dimensions of the polyhedra and the 


boundary conditions on the electron wave function will be 


49 





that the function be continuous between polyhedra and that 
it approach the zone boundary normally. On the average 
a monovalent metal will have only one electron per polyhedron. 
This is a result of both Coulomb repulSion and the effects 
of the Pauli exclusion principle. Thus each electron finds 
itself at a given instant confined to a atomic polyhedron 
mecluenced by the potentral oft the central 16n cone “Wagner 
and Seitz assumed a potential by Prokofjew [Ref. 19] to 
calculate the wave function and U,, the lowest state energy 
of this one elect.con system. 

Since the atomic polyhedron is not very different from 
a sphere, they assumed that it may be replaced by a sphere 
of equal volume. The radius of this sphere is Loe defined 
by = Tr =O where 2 is the atomic volume. 


The radialwwave functhen R = rl (r) so thatthe baund— 


ary condition on Y¥ which is =4 _. = 0 becomes 
Je ge 
aR _R (A2.2) 
SI alle) Os! Coa 


A convenient way to obtain the dependence of the energy 

on the radius of the sphere is to solve the radial Schrodinger 
equation for the wave function at arbitrary energies Uo: 

After obtaining the wave function, R(r), a tangent drawn 

from the origin to the curve will satisfy the boundary 
eemecition (A2~2)>°at two radi1l. “Thus Erom cach radial Eune@eron 
two points on the U(r) curve are determined. Calculation 
Serseveral Tadial funerron aliows the U(r) curve to pe 


mapped out. For sodium such a curve has a minimum and the 


50 








lattice would be stable without considering other contribu- 
tions to the energy. 

In fact however, the electrons in the metal are dis- 
tributed over states above Use the energy of the lowest 
State, by the Fermi distribution function. Thus to the 
energy of the lowest state must be added the average ad- 
ditional electron kinetic energy. This energy is deter- 


mined from the fermi energy 


2 
U, = = (22) AY 1, (ZB) 


where the average kinetic energy of an electron is 


ali 
Up aie Up (A2.4) 


These are the primary contributions to the cohesive energy 
of sodium. However, sodium has a relatively small ion 

core size so that the increase in the energy of the lattice 
from core overlap forces is not Significant. The larger ion 
core size in copper does not allow this assumption and the 
repulsive exchange energy between ions must also be con- 


Sidered. This energy is usually assumed to have the form 


_ sw 
U, = 5 (A2.5) 


where W is the repulsive bond energy of a single pair of 


ions and § is the lattice coordination number. W is usually 


taken to be of the form 


Wemper? © (AOmG) 
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Since op is a small fraction of the interatomic distance 
only first and second nearest neighbors are usually con- 
sidered to take part in the interaction. Other energies 
which are small and do not contribute substantially to 
the cohesive energy are the Van der Waals energy and a 
correction to the lowest state energy. This correction 

U. is the result of polyhedron-polyhedron interaction and 


L 


me che sum Of two terms. 
U. = U. - U (AZ. 7) 


where Up is the energy of a lattice of positive point 


charges immersed in a negative electron sea and U. in the 
energy of a positive point charge in a sphere of negative 
charge. 

The fact that a certain energy contributes signifi- 
cantly to the cohesive energy does not imply that it will 
Significantly affect the elastic constants. In obtaining 
the elastic shear constants of copper from the energies 
outlined above Fuch's used only the coulomb attraction 


energy of U, in equation (A2.7) and the overlap repulsion 


E 
energy of equations (A2.5) and (A2.6). 

The reasons the other energies do not contribute to 
the shear constants are as follows: 

a) For the equilibrium lattice spacing U, (r) is close 


to an inflection point and its second derivative 


1S close touzero. 


De 





b) For monovalent metals the Fermi surface is near 


the center of the Brillouin zone and is approxi- 


mately spherical. Thus shear strains do not alter 


the™average electron energy. 


c) The Fermi energy and the spherical contribution 


Oru and U are volume dependent energies. 


ra SY 


However, the strains chosen to define C and C' 
conserve the volume and thus these energies do 
not contribute to the to the elastic constants. 

ad) The Van der Waals energy makes only a minor 

GOnteE bDUutLIOon to Gand Cc’. 


Thus Fuch's retained only two contributions to the energy 


(A2.8) 


mimere U_ is the electrostatic 


E 
of point changes + e imbedded 


distribution of -e per atomic 


the exchange repulsion energy 


potential energy of a lattice 
in a uniform negative change 


volume and the energy U, is 


between closed shell ion 


cores. 


For a bcc lattice U,, 


20] to do the lattice sum and the electro- 


can be evaluated using the Ewald 
method [Ref. 
static contributions to the elastic constants are then 
calculated using equation (A2.1) so that 


2 2 


= = 2 = 
Slee OE) =r Cee Cee a 


(A2.9) 


where a is the lattice constant. 
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Similarily ionic stiffnesses Cy and C. are given by Fuch's 


I 


in terms of the energy W of equation (A2.6) 
Oe a ae SO tea A DO 
I 3 ya 2 i 2° 2 
Cc. = +n [r?2w'') + 2n [rw']) +N [rw'] (A210) 
9 1 1 9 : a | 2 ge 


Newe allowance has been made that the nearest neighbor in-= 
teraction a may differ from the next nearest neighbor 
interaction ite Primes indicate differentiation and the 
subscripts on the brackets denote evaluation at a distance 
r between ions corresponding to first or second nearest 
neighbors. 


The total shear constants are then 


= ina ' ' 
S Cr ata Ce. e Cr + Ce (A2.11) 


In practice Cr and Cr are usually determined by cal- 


culating CE and Ce and subtracting this from the measured 
constants. Featherston and Neighbours (Ref. 21] have done 
this for V, Nb, Ta, Mo, W and Fe. They note, in agreement 
with Isenberg, that for all cases except Fe the term 
- N(rwi), in Cy is negative and so must be offset by the 
second nearest neighbors. 

The application of this analysis to dilute alloys of 
copper was done by Neighbors and Smith. By introducing 
parameters to accomodate the effects of solute atoms on 


wre Clectron to atom ratio and the ionic overlap Contribue 


tion and by correcting for changes in the stiffnesses due 
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to changes in the lattice constant with solute they calculate 
Sorrections (toc ane. C = 

Further extensions of Fuch's theory to the transition 
elements and their alloys has encountered the difficulty 
of dealing with the status of the d-electrons. Such 
electrons are more tightly bound than s electrons, they 
thus occupy a relatively narrow band with a correspondingly 
higher density of states than s electrons. Further, the 
Bavave fumetions are more directional in character. ~The 
multi-valency of these metals and the changes in e/a with 
alloying will certainly cause changes in the fermi surface 
which could result in an electron gas contribution to the 
shear constants. The antiferromagnetic array that Zener 
and Isenberg predicted using Fuch's work has only been 


observed experimentally for chromium. 
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